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Abstract 

In  terms  of  fuel  cell  steady-state  performance  modelling,  many  electrical  models  have  been  developed  either  from  a  theoretical  point 
of  view  or  from  an  empirical  point  of  view.  The  model  described  in  this  article  is  from  the  empirical  point  of  view  approach.  This  model 
enables  to  simulate  both  fuel  cells  and  electrolysers  V-J  curves  (cell  voltage  versus  current  density)  in  typical  conditions.  This  model 
is  particularly  adapted  to  regenerative  fuel  cell  (RFC)  simulation.  It  is  a  four  degree-of-freedom  model  and  it  is  convergent  near  zero 
current.  It  depends  on  the  stack  temperature  and  the  oxygen  partial  pressure.  The  regions  where  mass  transfer  limitations  occur  have  not 
been  modelled,  because  they  are  usually  avoided  for  efficiency  or  thermal  reasons.  The  parameters  have  been  fitted  with  a  4kWe  proton 
exchange  membrane  fuel  cell  (PEMFC)  and  a  3.6  kWe  electrolyser.  The  electrical  equations  and  the  experimental  data  are  well  correlated. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  Centre  d’Energetique  of  Ecole  des  Mines  de  Paris 
at  Sophia  Antipolis  has  gained  experience  in  fuel  cell  (FC) 
studies  for  10  years,  and  particularly  in  proton  exchange 
membrane  fuel  cell  (PEMFC).  Since  the  development  of 
the  first  fuel  cell  stack  test  bench,  many  works  have  been 
achieved  in  terms  of  fuel  cell  integration,  hydrogen  genera¬ 
tion  or  storage  and  fuel  cell  modelling.  For  testing  purposes, 
a  dedicated  fuel  cell  stack  test  bench  is  available,  capable  of 
testing  fuel  cells  up  to  15  kWe  in  various  system  configura¬ 
tions. 

The  laboratory  is  also  doing  research  concerning  hybrid 
systems  made  up  of  a  renewable  energy  generator  (photo¬ 
voltaic  field  or  PV),  a  back-up  unit  (diesel  engine)  and  a 
storage  system  (batteries,  fuel  cell/gas  storage/electrolyser). 

For  these  projects  two  test  facilities  have  been  set  up.  The 
first  one  is  a  PV-battery-diesel  test  bench  and  the  second 
one  is  a  PV-electrolyser-fuel  cell  test  bench. 

This  latter,  constructed  within  the  frame  of  the  PVFC-SYS 
European  project  (ERK5-CT1999-00017)  is  an  autonomous 
electricity  generator  including  a  PV  field  (3.6  kWp),  an  elec¬ 
trolyser  (3.6  kW)  splitting  water  into  hydrogen  and  oxygen, 
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a  gas  storage  unit  (4Nm3  of  H2  and  2Nm3  of  O2)  and  a 
PEMFC  (4  kWe)  to  generate  electricity  during  low  sun  shine 
periods. 

This  test  bench,  described  in  different  articles  [1-3]  en¬ 
ables  us  to  validate  the  model  of  each  component  and  finally 
of  the  complete  system. 

Different  approaches  exist  to  simulate  the  electrical  char¬ 
acteristics  of  a  FC.  There  are  roughly  two  kinds  of  fuel  cell 
models. 

The  first  one  is  the  theoretical  model  (also  named  mecha¬ 
nistic  model),  which  describes  the  electrochemical  reactions 
occurring  in  the  cells.  In  the  model  of  Metkemeijer  [4],  the 
cell  voltage  is  the  thermodynamical  voltage,  minus  the  dif¬ 
ferent  overvoltages  due  to  the  ionic  transfers  at  the  anode  and 
at  the  cathode,  the  resistive  losses  and  the  material  transfer 
at  high  current  density.  Each  term  of  the  equation  depends 
on  the  stack  temperature  and  the  partial  pressure  of  hydro¬ 
gen  and  oxygen.  This  approach  requires  the  knowledge  of 
nine  parameters,  which  are  difficult  to  determine. 

The  model  of  Amphlett  [5],  which  is  largely  quoted  in 
the  literature,  is  also  based  on  Nernst  and  Tafel  equations. 
It  considers  all  physical  parameters  in  the  system  (effective 
pressure  of  oxygen  and  hydrogen,  temperature,  concentra¬ 
tion  of  oxygen,  hydrogen,  water,  proton).  Since  all  these 
parameters  cannot  be  identified,  the  authors  use  empirical 
means  to  estimate  their  values  [6] . 
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Since  then,  it  has  been  improved  and  now  allows  to  take 
into  account  the  fuel  cell  parameters  (such  as  active  area 
and  membrane  thickness)  and  the  voltage  degradation  with 
membrane  ageing  [7,8].  This  model  has  to  consider  empir¬ 
ical  and  semi  empirical  parameters,  like  ohmic  overvoltage 
and  the  equilibrium  water  content  of  the  membrane. 

In  the  second  method,  a  semi  empirical  equation  is  used, 
of  which  parameters  are  calculated  through  an  identification 
process  with  experimental  data.  The  terms  of  this  equation 
come  from  Nernst,  Tafel  and  Ohm  laws.  This  kind  of  model 
does  not  pretend  to  be  universal  but  it  is  much  simpler  than 
theoretical  ones. 

The  model  of  Chamberlin  and  Kim  [9]  describes  the  cell 
voltage  depending  on  the  current  density,  with  five  parame¬ 
ters  (Eq.  (1)). 

Vcell(T)  =  Eo  —  b\n(J)  —  R  x  J  —  m  exp (nJ)  (1) 

where  Vceii  is  the  cell  voltage  (V),  J  the  current  density 
(A cm-2),  Eo  the  open  circuit  voltage  (V),  b  the  Tafel’s 
parameter  for  the  oxygen  reduction  (V),  R  the  resistance 
(£2  cm2),  m,  n  is  the  diffusion’s  parameters  (V  and  A-1). 

Each  term  of  Eq.  (1)  is  dominant  in  each  region  of  the 
V-J  characteristic  (Fig.  1).  In  region  1,  the  voltage  decreases 
drastically  due  to  the  oxygen  electrochemical  activation  re¬ 
actions  (logarithm  term  has  the  main  influence).  In  region 
2,  the  curve  is  roughly  linear  (resistive  losses).  The  region 
3  corresponds  to  the  diffusion  losses  (exponential  term). 

In  Kim’s  model,  the  five  parameters  (Eo,  b ,  R ,  m,  n)  de¬ 
pend  on  the  temperature  and  pressure  in  the  cell  and  the  oxy¬ 
gen  partial  pressure.  They  also  depend  on  the  stack  itself, 
so  that  it  cannot  be  transposed  to  another  fuel  cell  without 
a  new  identification  of  the  parameters. 

This  model  fits  very  well  with  experimental  results. 
Hamelin  has  compared  the  Amphlett’s  model  and  the  Kim’s 
one  [12].  It  appears  that  both  predicted  voltages  are  very 
close  to  each  other  and  that  they  both  differ  slightly  from 
experimental  cell  voltage  in  the  same  regions  (low  and  high 
current  densities). 


However,  it  is  a  real  issue  to  simulate  the  fuel  cell  be¬ 
haviour  in  the  first  zone  (i.e.  for  low  current  densities).  In¬ 
deed,  the  Kim’s  model  is  divergent  close  to  J  =  0. 

This  problem  particularly  appears  when  modelling  a  re¬ 
generative  fuel  cell,  a  device  that  can  be  used  either  as  an 
electrolyser  or  as  a  fuel  cell,  because  equation  (1)  diverges 
each  time  the  device  passes  continuously  from  one  opera¬ 
tion  to  the  other  one. 

To  solve  this  problem,  we  propose  a  new  model  passing 
through  the  divergence.  Described  in  chapter  II,  the  model 
is  validated  in  chapter  III  for  a  fuel  cell  and  an  electrolyser 
running  in  regions  1  and  2. 

To  simulate  an  electrolyser,  the  same  electrical  equation 
is  used,  changing  the  value  of  the  parameters. 

The  model  is  hence  able  to  simulate  a  regenerative  fuel 
cell  throughout  its  complete  behaviour. 


2.  Description  of  the  model 

The  proposed  model  is  very  close  to  the  classical  semi 
empirical  model  of  Kim  with  the  advantage  to  be  conver¬ 
gent  for  all  current  density.  Furthermore,  it  is  a  purely  em¬ 
pirical  model.  The  electrical  equation  as  well  as  the  matlab 
program  have  been  developed  in  the  laboratory  in  1999.  It 
is  a  four-degree-of-freedom  model. 

This  model  does  not  take  into  account  the  diffusion  area 
(region  3  in  Fig.  1).  Like  in  Kim’s  model,  it  is  possible  to 
add  an  exponential  term  to  tackle  this  issue.  The  authors  did 
not  want  to  damage  the  fuel  cell  by  using  it  in  this  region. 
Fuel  cells  are  usually  used  in  a  region  where  no  thermal  or 
diffusion  problems  occur  for  all  temperatures  (normally  for 
cell  voltages  >  0.5-0. 6  V). 

2.7.  Electrical  equation 

To  solve  the  convergence  problem,  we  propose  to  change 
the  logarithm  term  into  (1/ln  J),  close  to  0  for  J  =  0.  The 
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following  equation  is  obtained: 

Veen  (J)  =  E+-^—-cJ  (2) 

In  (dJ) 

where  Vceii  is  the  cell  voltage  (V),  /  the  current  density 
(A  cm-2),  E ,  b ,  c,  d  the  parameters  (V,  V,  Q  cm-2,  cm2  A-1). 
The  other  hypotheses  are  summed  up  in  Fig.  2. 
Parameter  E  has  a  physical  meaning:  it  is  the  open  circuit 
voltage. 

The  curve  of  the  cell  potential  (V)  versus  the  current  den¬ 
sity  (J)  has  an  inflexion  point  for  J  —  Jd  and  at  this  point, 
the  slope  of  the  tangent  is  noted  —  A;  it  is  the  sum  of  acti¬ 
vation  and  resistive  losses  (see  Fig.  1). 

From  the  two  equations:  V"CQ\i(Jd)  =  0  and  V^el l(Jd)  = 
—  A,  we  obtain  two  new  equations: 


and 


b 


Replacing  this  terms  in  Eq.  (2),  we  finally  end  up  with: 

y“"<J)  -  £  +  2  +  (i  - A)  *  '  <3) 


where  E  is  the  open  circuit  voltage  (for  J  =  0),  /j,  A  the 

other  parameters  of  the  model  (V,  A  cm-2,  £2  cm2). 


2.2.  Four  parameters  to  be  determined 


The  four  parameters  ( E ,  /j,  b  and  A)  have  to  be  deter¬ 
mined,  fitting  the  modelling  curve  to  the  experimental  data. 
These  parameters  depend  on  the  cell  temperature  and  on 
the  oxygen  partial  pressure.  The  influence  of  the  hydrogen 
partial  pressure  can  be  neglected  as  long  as  the  hydrogen 
fraction  on  the  anodic  FC  side  is  superior  to  20%  [10]. 

Due  to  the  temperature  and  pressure  dependence  of  the 
thermo  chemical  potential  and  the  Butler  Volmer  equation, 
the  influence  of  T  and  po2  on  the  four  parameters  has  the 
following  form:  K\  +  K2T  +  K^T  ln(po2). 

Therefore,  for  each  parameter,  three  constants  have  to 
be  determined  (Eq.  (4))  in  order  to  describe  the  T  and  P 
dependence. 


"  E 

Jd 

b 

_  A 

T  in  K,  po2  in  bar. 

To  do  so,  a  minimum  of  16  data  is  required:  four  couples 
(/,  Veen)  for  four  pairs  { po2 ,  T}  where  po2  and  T  must  be 
at  least  at  two  different  levels. 

In  case  of  the  electrolyser,  the  oxygen  partial  pressure  cor¬ 
responds  to  the  operating  pressure  of  the  electrolyser  (pure 
oxygen  is  produced). 

In  case  of  the  fuel  cell,  when  pure  oxygen  is  used,  the  par¬ 
tial  pressure  corresponds  to  the  system  pressure.  When  am¬ 
bient  air  is  used,  the  partial  pressure  of  the  oxygen  depends 
on  the  oxygen  molar  fraction  and  on  the  oxygen  utilisation 
(inverse  of  stoichiometric  factor). 

2.3.  Matlab®  model 

To  work  out  the  four  parameters  ( E ,  Jd ,  b  and  A),  the 
mathematical  method  used  is  the  least  square  method  for 
non-linear  problems.  Knowing  the  [4x3]  matrix  of  param¬ 
eters,  the  cell  voltage  is  a  function  of  the  current  density,  the 
temperature  and  the  oxygen  partial  pressure. 


E  i 

E2 

E  3 

Jd  1 

Jdl 

Jd3 

b\ 

b2 

h 

^2 

2^3 

1 

T 

_T\n(po2)  _ 

(4) 


3.  Experimental:  fitting  of  empirical  equation 

The  fuel  cell  and  the  electrolyser  were  tested  on  the  last 
test  bench  described  in  the  introduction. 

The  electrolyser  is  a  3.6 kW  alkaline  one  from  Hydrogen 
Systems,  which  has  16  cells  of  300  cm2  connected  in  series. 
The  periphery  has  been  developed  and  adapted  to  the  power 
and  to  the  application.  The  cell  temperature  is  calculated 
as  the  average  between  the  two  inputs  and  two  outputs  of 
the  electrolyte  in  the  cells.  The  pressure  is  measured  in  the 
hydraulic  periphery. 
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Fig.  3.  Experimental  data  and  simulation  curves  of  the  electrolyser — influence  of  the  temperature  (P  =  9.6  bar). 


The  FC  system  (described  in  [11])  has  been  optimised 
to  our  specific  application.  The  proton  exchange  membrane 
stack  from  DeNora,  is  made  up  of  18  cells  of  900  cm2. 

The  dedicated  FC  system  enables  to  operate  for  variable 
oxygen  (and  hydrogen)  partial  pressures  and  at  variable  tem¬ 
peratures. 

The  system  pressure  can  vary  from  1  to  4  bar  (abs).  To 
obtain  the  variable  oxygen  and  hydrogen  partial  pressures, 


oxygen  and  hydrogen  were  diluted  with  nitrogen.  In  order 
to  avoid  mismeasurements  due  to  nitrogen  migration  across 
the  membrane  [10],  hydrogen  and  oxygen  partial  pressures 
were  always  balanced. 

In  the  FC,  the  cell  temperature  is  not  available.  Usu¬ 
ally,  the  temperature  of  the  air  output  is  selected  because  it 
represents  most  precisely  the  electrode-membrane-assembly 
temperature.  In  our  approach,  we  try  to  simulate  empir- 


Fig.  4.  Experimental  data  and  simulation  curves  of  the  electrolyser — influence  of  the  pressure  ( T  =  52-56  °C). 
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Fig.  5.  Experimental  data  and  simulation  curves  of  the  FC — influence  of  the  temperature  ( po2  =  1.4  bar). 


ically  the  complete  FC  system  and  then  the  temperature 
selection  is  not  a  crucial  point.  For  technical  reasons,  the 
temperature  used  in  the  electrical  equation  is  the  tempera¬ 
ture  of  the  inlet  cooling  water  (with  an  outlet  cooling  wa¬ 
ter  temperature  maximum  5  °C  above  inlet  temperature). 
For  this  stack,  the  optimal  temperature  has  been  evaluated 
to  35  °C. 


3.1.  Results 

Each  constant  has  been  determined  for  the  electrolyser 
and  for  the  FC  as  a  function  of  temperature  and  oxygen 
partial  pressure  by  using  the  experimental  results  obtained 
with  the  FC  and  the  electrolyser  previously  mentioned  (see 
Figs.  3-6).  Although  the  experimental  data  are  a  bit  dis- 


Fig.  6.  Experimental  data  and  simulation  curves  of  the  FC — influence  of  the  pressure  (T  =  35  °C). 
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Table  1 


Constants  of  the  electrical  model:  set  of  values  for  the  electrolyser  (left)  and  for  the  fuel  cell  (right) 


Electrolyser 

Fuel  cell 

Jd 

0.355 

—  1.078E— 3 

5.232E— 4 

-0.868 

—  1.762E— 3 

4.9E-4 

E 

1.025 

1.24E-3 

-1.05E-4 

0.796 

4.661E— 3 

-3.22E-4 

b 

-6.125 

1.57E-2 

— 3.245E— 4 

-1.926 

2.616E— 2 

—  1.489E— 3 

A 

-1.819 

4.83E-3 

—  1.248E— 4 

3.892 

— 7.447E— 3 

— 2.696E— 4 

Jd'.  Am  2;  E:  V;  b :  V;  A:  Qm2. 

persed  due  to  the  noise  in  the  measurement,  the  fitting  pro¬ 
tocol  enables  to  select  the  most  representative  curve.  The 
identification  process  leads  to  the  two  following  [4x3]  ma¬ 
trixes  (Table  1). 

3.2.  Discussion 

To  validate  the  model,  we  can  observe  the  influence  of  the 
temperature  and  the  oxygen  partial  pressure.  We  will  point 
out  the  limits  of  the  model  validity,  but  before  it  should  be 
born  in  mind  that  the  model  prediction  is  only  valid  for  the 
specific  equipment  and  inside  the  temperature  and  pressure 
range  used  to  fit  the  parameters  (usually  real  operating  lim¬ 
its).  For  the  electrolyser,  the  experimental  pressure  is  be¬ 
tween  1.7  and  9.9  bar  (abs)  and  the  temperature  varies  be¬ 
tween  31  and  63  °C.  For  the  fuel  cell,  the  oxygen  and  hy¬ 
drogen  partial  pressures  range  from  0.7  to  1.5  bar.  The  tem¬ 
perature  range  is  15-36  °C. 

3.2.1.  Temperature  influence 

The  Figs.  3  and  5  plot  the  V-J  characteristics  of  the  elec¬ 
trolyser  and  the  FC  at  fixed  pressure  for  different  tempera¬ 
tures.  As  expected,  the  temperature  increase  induces  a  better 
performance  of  both  the  electrolyser  and  the  FC  in  terms  of 
voltage  efficiency. 

3.2.2.  Pressure  influence 

The  pressure  has  an  important  effect  on  the  FC  perfor¬ 
mance  (Fig.  6)  but  not  on  the  electrolyser’ s  one  (Fig.  4). 
More  precisely,  the  high  pressure  has  an  influence  only 
on  the  purity  of  gases  produced  by  electrolysis,  but  not 
on  the  reaction  efficiency.  A  model  of  the  gas  purity  will 
be  added  later,  considering  the  flux  crossing  through  the 
membrane  depending  on  the  electrolyser  pressure.  As 
predictable,  the  oxygen  partial  pressure  increase  induces 
a  better  performance  of  the  FC.  As  the  system  pres¬ 
sure  is  3  bar  absolute,  the  oxygen  concentration  varies 
between  25%  (0.7  bar)  and  50%  (1.5  bar).  For  safety  rea¬ 
sons,  we  did  not  measure  in  regions  higher  than  50%  of 
oxygen. 

3.3.  Conclusions  concerning  the  electrical  model 

The  proposed  model  is  a  first  step  model.  It  is  a  satis¬ 
factory  model  to  simulate  an  electrolyser  or  a  FC  for  low 
and  medium  current  density.  This  model  is  innovative  due  to 


its  convergence  at  zero  current  density  and  its  mathematical 
coherence. 

The  model  shows  a  good  fit  of  the  empirical  equation  to 
the  experimental  data  and  the  influence  of  the  pressure  and 
temperature  parameters  is  correct  compared  to  the  physical 
behaviour  of  the  two  electrochemical  components. 

4.  System  modelling 

The  electrical  model  is  only  a  part  of  the  complete  elec¬ 
trolyser  or  fuel  cell  model.  As  seen  above,  the  electrical 
model  needs  two  inputs:  the  temperature  and  the  oxygen  par¬ 
tial  pressure.  Two  sub  models  are  essential  to  entirely  simu¬ 
late  the  electrochemical  devices.  The  first  one  is  the  thermal 
model,  allowing  to  simulate  the  evolution  of  the  tempera¬ 
ture  during  the  component  operation.  The  second  one  is  the 
gas  management  model  calculating  particularly  the  oxygen 
partial  pressure. 

The  Fig.  7  shows  the  electrolyser  model,  programmed 
with  Matlab/Simulink®,  made  of  three  main  modules  allow¬ 
ing  to  calculate  the  electrical  operating  points,  the  temper¬ 
ature  and  the  pressure  of  the  electrolyser  system.  The  fuel 
cell  system  uses  the  same  program  architecture.  The  elec¬ 
trolyser  model  has  two  inputs  (power  and  ambient  tempera¬ 
ture)  whereas  the  FC  model  has  three  inputs  (power,  ambient 
temperature  and  oxygen  partial  pressure). 

Other  sub  modules  could  be  of  great  importance  to  com¬ 
plete  the  simulation  depending  on  its  application.  Indeed,  the 
simulation  of  the  intrinsic  consumption  of  the  peripheries  is 
required  to  perform  the  evaluation  of  the  complete  system 
efficiency,  which  can  be  done  by  a  dedicated  sub  module. 
In  addition,  in  some  particular  applications,  transient  phe¬ 
nomena  could  occur,  which  can  be  managed  by  another  sub 
module. 

4.1.  Thermal  model 

To  estimate  the  temperature  evolution  during  operation, 
the  following  equation  can  be  used. 

fd0\ 

Cp  I  —  I  —  Pth  —  ^ext  —  ^ech  —  ^gaz 

0=T-Ta 

Pth  =  3znc(U  —  Utn)I  (thermal  dissipative  power) 
cpQXt  =  h0  (heat  flow  lost  in  the  atmosphere) 
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Fig.  7.  The  complete  electrolyser  model  (without  transient  module)  with  Matlab/Simulink®- 


<Pg az  =  ±(Cp(H2)F(H2)  +  CP(02)F(02))  x  # 

(heat  flow  lost  in  gas)  (+for  the  electrolyser, 

—  for  the  fuel  cell) 

where  Cp  is  the  thermal  capacity  of  the  device  (JK-1),  T 
the  temperature  (K),  Ta  the  ambient  temperature  (K),  nc  the 
number  of  cells,  U  the  cell  voltage  (V),  Utn  the  thermoneu¬ 
tral  cell  voltage  (=1.48  V)  based  on  the  higher  heating 
value  of  hydrogen,  I  the  current  (A),  h  the  heat  transfer 
coefficient  (WK-1),  (pech'  energy  flow  lost  in  a  heat  ex¬ 
changer  (WK-1),  Cp(i)  the  molar  heat  capacity  of  the  gas 
i  (Jmol-1  K-1),  F(t)  the  molar  flow  of  the  gas  i  (mols-1). 

The  thermal  dissipative  power  Pth  is  calculated  knowing 
the  electrical  operating  point  (current-voltage).  The  two 
coefficients,  Cp  and  h ,  are  fitted  with  experimental  data. 
Depending  on  the  system  component,  the  heat  transfer  coef¬ 
ficient  h  can  depend  on  the  current  (usually  for  the  alkaline 
electrolyser  due  to  the  electrolyte  circulation  proportional 
to  the  current). 

This  thermal  module  can  contain  a  heat  exchanger  (<pech) 
model  allowing  to  regulate  the  device  temperature  and  to 
calculate  the  heat  production  of  the  system. 

4.2.  Pressure  model 

The  gas  management  module  enables  to  calculate  gas 
production  or  consumption  and  the  periphery  pressure  level 


during  operation.  The  gas  purity  of  the  electrolyser  products 
will  be  simulated  shortly. 

To  obtain  the  gas  production  or  consumption,  the  two  next 
equations  are  used: 


ncI 

^gas  —  —  x  Of( EL) 

nF 


ncI 

nF 


1 

x  — 
OF 


(FC) 


where  Pgas  is  the  gas  consumption  for  the  fuel  cell  and  gas 
production  for  the  electrolyser  (mols-1),  nc  the  number  of 
cells,  I  the  component  current  (A),  ijf  the  faraday  efficiency 
(%),  F  the  faraday  constant  (96485  C  mol-1),  n  the  number 
of  moles  of  electrons  transferred  per  mole  of  water  (n  =  2 
for  hydrogen,  n  =  4  for  oxygen). 

The  pressure  calculation  of  the  electrolyser  can  be  in¬ 
teresting  according  to  the  electrolyser  system  and  its  ap¬ 
plication,  although  it  has  few  influence  on  its  electrical 
behaviour.  During  the  start  up  of  the  device,  the  periph¬ 
ery  pressure  increases  slightly  until  the  maximal  operat¬ 
ing  pressure.  For  our  system  (operating  pressure  at  lObar), 
it  usually  lasts  about  20  min  but  it  could  be  longer  for  a 
high-pressure  electrolyser.  Moreover,  pressure  has  an  influ¬ 
ence  on  the  gas  purity.  Especially,  hydrogen  tends  to  ac¬ 
cumulate  in  the  oxygen  circuit,  which  may  lead  to  explo¬ 
sive  mixtures.  Electrolyser  pressure  calculation  is  therefore 
essential. 

As  for  the  fuel  cell  system,  the  oxygen  partial  pressure  is 
fixed  at  the  start  up  of  the  device  and  remains  constant  during 
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operation,  the  gas  management  module  calculates  only  gas 
consumption. 

To  calculate  the  gas  management  unit  results,  some  pa¬ 
rameters  have  to  be  determined  experimentally:  faraday  ef¬ 
ficiency,  gas  flux  depending  on  the  periphery  pressure,  etc. 

5.  General  conclusions 

The  model  presented  here  is  a  novel  empirical  model  to 
precisely  calculate  the  major  part  of  the  V-J  characteristic 
of  a  fuel  cell,  an  electrolyser  or  a  reversible  fuel  cell.  With  a 
set  of  experimental  data,  the  parametric  electrical  equation 
can  be  easily  fitted  with  any  kind  of  fuel  cell  or  electrol¬ 
yser.  This  model,  which  is  mathematically  consistent  and 
convergent  near  zero,  allows  the  operator  to  interpolate  ex¬ 
perimental  results  of  his  studied  device  running  in  regions 
where  current  density  does  not  cross  the  mass  transfer  limi¬ 
tations.  This  model  will  be  soon  used  to  study  the  influence 
of  the  hydrogen  partial  pressure  on  the  FC  performance.  In 
case  of  validation,  the  model  would  enable  to  simulate  FC 
behaviour  supplied  with  different  fuel  inputs. 

The  electrical  model  is  an  essential  module  of  an  electro¬ 
chemical  device  simulation  but  many  other  points  are  im¬ 
portant  during  the  operation.  Each  module  proposed  in  the 
fourth  part  can  be  useful  depending  on  the  use  of  the  elec¬ 
trolyser  or  the  fuel  cell.  Nevertheless,  in  our  application  all 
these  models  are  needed  to  perfectly  simulate  the  complete 
system. 
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